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Radiowave propagation in the extremely low frequency (ELF) range of 3 Hz to 3 kHz has been 
successfully treated in the past by a type of waveguide mode theory. This theory simplifies at large 
distances where only the zero order mode need be considered. It is shown in this paper that for 
frequencies in the band of 5 Hz to 1 kHz the zero order mode predominates at any distance in excess 
of roughly 500 km. 

In a lower part of the ELF band (5 Hz to 100 Hz) more than one approximating form must be used. 
The commonly used asymptotic expansion applies only to a middle range of distances between source 
and receiver. When the receiver is within roughly one-sixth wavelength of either the source or its 
antipode then different approximations must be employed. These extensions of the theory are used 
in this paper to derive modifications of previous results for the terrestrial ELF noise power spectrum. 



1. Introduction 

The theory of extremely low frequency (ELF, 3 Hz to 3 kHz) radiowave propagation, developed 
by Schumann [1957] and Wait [1960a, b], predicts that the spherical cavity between the earth and 
the ionosphere will resonate at various frequencies for which the phase retardation in one circuit 
around the earth equals an integer number of complete cycles. This resonance phenomenon was 
observed experimentally in the natural noise background by Balser and Wagner [1960] over the 
"resonant" range of frequencies (5 to 50 Hz). These observations were fitted in some detail by 
Raemer [1961] who derived an expression for the power spectrum of the lightning flashes thought 
to be responsible. Improved explanations were given by Galejs [1961a, 1962] who had pre- 
viously developed a theory of ELF propagation (5 Hz to 1 kHz) for a more realistic model of ex- 
ponentially increasing conductivity with height in the isotropic ionosphere [Galejs, 1961b]. The 
earth's magnetic field was mostly ignored in these references and this paper will follow suit. The 
effect of such a static magnetic field has already been shown by Wait [1960b] to be largely 
negligible in the ELF band. 

Some additions to the existing theory are needed in the resonance range when the noise source 
is close to the observer or his antipode. The work of both Raemer and Galejs was compared to 
observations taken by Balser and Wagner [1960, 1962] in the vicinity of Boston, Mass. The 
lightning stroke noise sources in these cases were mostly far distant from the receiver or its 
antipode. Hence, in this sense, the results of Raemer and Galejs are correct as they stand. The 
point of this paper is that they cannot be translated directly to measurements performed at sites 
close to the great tropical thunderstorm areas such as the Amazon or the Congo River basins (or 
their antipodes). All of the work presented here applies to the frequency range of 5 Hz to 1 kHz 
in agreement with the work of Galejs [1961b]. 

The waveguide mode theory of radio propagation is useful at ELF [Wait, 1960b] where the 
earth and the ionosphere form the lower and upper boundaries, respectively. The distance along 
the earth's surface is measured by the angle subtended at its center. The individual modes in 
such a spherical system have a dependence on distance given by P v {n){~ cos 0), the Legendre func- 
tion of argument (—cos 0) and of complex degree v(n). In this paper the zero order mode (n = 0) 
is defined to be the mode of least attenuation in the ELF band. It is commonly stated that at large 
distances the zero order mode will predominate and the others can be ignored. This paper shows 
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that the critical distance at which this first occurs is roughly equal to 500 km; a value large com- 
pared to the height of the ionosphere and, to the first order, constant with frequency in the range 
of 5 Hz to 1 kHz. The mode cutoff point around 3 kHz is not treated. 

This last statement may seem to be at variance with previous ideas that one-sixth wavelength 
is a minimum distance. What Wait [1960b] has shown is that the normalized magnitude of the 
electric field varies in a linear manner with distance when the distance exceeds one-sixth wave- 
length and the frequency is in the range of 50 to 1600 Hz. However, the fact is that one-sixth 
wavelength also arises very naturally as the minimum distance needed to use a certain simplifying 
asymptotic expansion for P„(o)(— cos 6). For frequencies below 100 Hz there is a range of distances 
running from 500 km out to one-sixth wavelength where this first asymptotic form does not hold. 
But the zero order mode still predominates. Fortunately, a second expansion is available that does 
hold. Similarly, a third approximate form must be used when the radio wave source is within 
one-sixth wavelength of the observers' antipode. A second result of this paper is to present, 
under these latter two conditions, the necessary modifications to the formula of Galejs [1961a] 
for the terrestrial noise power spectrum. 

2. Background Results 

In support of the main points of this paper, it is useful to set down some previously developed 
approximations for the Legendre function P u ( — cos 6) and the associated Legendre function of, 
first order PJ(— cos 6). An approximation for the complex degree v{n) as a function of Ji(=l, 2, . . .) 
is also given. These are results that are needed in various parts of this paper. 

2.1. Approximate Forms for P v (—cos 0) and Pv 1 (—cos 6) 

The various mathematical approximations quoted here can be found either in Bateman 
Manuscript Project [1953] or Magnus and Oberhettinger [1949]. They are grouped according 
to the relative position of the noise source and the receiver. Further approximations are shown 
that apply more particularly for the zero order mode since the complex degree ^(0) is largely real 
[Galejs, 1961b]. 

The asymptotic expansions conventionally used in ELF propagation theory apply when the 
receiver is not too close to the source or its antipode. More explicitly, 



^-^ 9) ~ W^k Vdb C ° S K" + a/2>*r-«)-<ir/4)] 



P ' ( - COS 9) ' W^k) Vdb C ° S K" + <1/2)X-- » + <T/4)] 

if\v\6> 1 and \v\{tt-0)>\ and (1) 

where Y(z) is the gamma function of z. 

Each of these is the first term of a trigonometric expansion. If the Sterling formula is used to 
approximate the gamma functions, the asymptotic expansions reduce to 



P v ( - cos 0) - V2/(*; + (l/2))7rsin<9 cos [{v + (1/2))(tt - 0) - (tt/4)] 



PJ( - cos 9) - V2(i; + (l/2))/7r sin 6 cos [(v + (1/2))(tt - 6) + (tt/4)] 

if\v\6>\ and \v\ (tt-0)>1. (2) 
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These last approximations are most appropriate in the case of the zero order mode where v(0) 
is largely real. 

When the source is close to the antipode of the receiver an expansion of P v {— cos 0) for small 
{tt — 0) is required. Probably the most rapidly convergent is a series, due to MacDonald, in terms 
of Bessel functions [Bateman Manuscript Project, 1953]. However, for the computations in 
this paper a better choice is based on the hypergeometric series definition of the Legendre function, 

P„(-cos0)^ 1 

PJ(- cos 6) - - v(y + 1) cos ((9/2) 

if\v\(7T-d)<\. ( 3 ) 

Finally, when the sources are near the receiver itself a third approximation must be used to properly 
account for the singularity there. Materially different appearing expressions apply, depending 
on whether the order is zero or positve [Bateman Manuscript Project, 1953]: 

PA- cos 0) - (sin (vn)lir) [2 log (sin (0/2)) + C + 2\jj(v + 1) + n cot (ptt)] 

Pl(- cos 6) - - (sin (v7r)lir)lsin (0/2) 

if\v\6< 1 and (4) 

where C = Euler's constant = 0.5772 — 
\\f{z) = logarithmic derivative of the Gamma function 

= d(\ogV(z))ldz = -C + ^[(ll(n+l))-(ll(n + z))]. 

tt=0 

These last results complete the approximations needed to cover the entire range of polar angle. 

2.2. Higher Order Roots of the Modal Equation 

Wait [1960a] has already shown in the ELF band for a homogeneous ionospheric layer that 
only the zero order mode is freely propagating. The higher order modes (n ^ 1) are below cutoff 
and thus heavily attenuated with distance. The same result can be obtained, using the zero order 
mode solution by Galejs [1961b], for a more realistic, exponentially increasing conductivity of 
the ionosphere. (The exponential model had been previously treated by Wait [1960a, b]. The 
solution, due to Galejs, represents an extension to that work.) 

The complex degree v(n) is related to a mode factor S n by the relation 

v{v+l) = {kaSn? (5) 

where 

/c = free space, plane wave propagation constant 
a = radius of the earth. 

The so-called modal equation is written in terms of S n to satisfy boundary conditions at the earth 
and at the ionosphere. In the ELF band Galejs showed for the exponential ionosphere that So is 
largely real and close to unit magnitude, i.e., Re(S ) > 10 Im(S ) and S — 1 for frequencies from 
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5 Hz to 1 kHz. He observed that this solution can be used to define an equivalent homogeneous 
ionosphere that will have an identical surface impedance Z\. 

It can be said that S n is a function of Z* which, in turn, is related to S which is already known. 
Using the basic relations of Wait the end result is 



S I = 1 - (irnfflK? [ 1 + Vl+(2kh/Trn) 2 (l-Sl)] 2 (6) 

where h = height of the ionosphere above the earth. Detailed examination of numerical values for 
So [Galejs, 1961b, 1962] shows that the square root can be expanded with little error. The param- 
eter (2khlirn) 2 {l—Sl) has a maximum value of 0.1 at the upper frequency limit of 1 kHz. 
After performing the expansion the higher order factors are approximated by 

S 2 --(™/£/i) 2 + (2S 2 -l). (7) 

Solving for the higher order complex degree (5) shows that to a close approximation 

v{n) - -jiiranlh) +j{k 2 ahl2>nn) (2S 2 - 1) - -jinan/h). (8) 

In this form, the first constant term on the right predominates (as shown) for frequencies from 
5 Hz to 1 kHz. The higher order modes (n ^ 1) have a complex degree that is largely imaginary. 
Thus, it follows that they are rapidly attenuated with distance. 



3. Critical Distances 

The previous results are used in this section to help prove that the zero order mode begins to 
predominate at distances in excess of 500 km. Also, the ranges of polar angle are specified where 
various approximate forms for this zero order mode apply. 

3.1. Minimum Distance for Predominance of the Zero Order Mode 

The vertical electric field at the earth's surface for a vertical current element at the surface is 
given by [Wait, 1960a,b], 



E r = (r)Ids/2kha 2 ) 



v(0) (i/(0) + l)/\<o)( ~ cos 0)12 sin v(0)tt 
+ ]T i/(i/+1)P„(-cos 0)/sin vir 



(9) 



where 6 = polar angle subtended by the source and the observer at the center of the earth, 
17 = characteristic impedance of free space — 120 7rfl, 
k = free space, plane wave propagation constant, 
Ids = moment of the vertical element, 
a = radius of the earth, 
h = height of the ionosphere above the earth. 

Near the source all of the modes are required. At larger distances only the zero order mode has 
any significance since all the higher order modes are below cutoff. 

The cutoff effect can be clearly seen by the following approximation which applies at larger 
distances for the higher order mode term in the series for electric field. 



|i<i/ + l)P,(-cos 0)/sin vtt\ - (rran\K?^\/2\TT sin 6 exp (- (7ranlh)6) (10) 
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for n ^ 1 



if\v\0>l and \p\(tt-6)> 1. 

(This result uses those of (2), (5), (7), and (8).) The parameter (tt an/h) has a minimum value of 
about 200 for an ionosphere of maximum height, h = 100 km. Hence, the exponential reduces the 
amplitude to negligible proportions. 

The same type of approximation for the zero order mode term is 



|i^+l)/V-cos 6)12 sin vtt\ - v(y + 1) cos[(i/ + (1/2))(tt- 6)- (7r/4)]/(sin i>tt)V(2i/ + 1)tt sin 6/ 

(ID 



for n = 



i/|i/|0>l and H(tt-0)> 1. 



The degree ^(0) has close to unit value. 

By comparison to the higher order mode (see (10)) the zero order mode predominates at any 
large distance. Even its nulls are blunted by the slight propagation losses of the zero order mode. 

As the distance (polar angle) is reduced, the logarithmic increase of P„(— cos 6) (see (4)) even- 
tually comes into play where all modes become important. In order for the zero order mode to 
predominate it is required that \v{n)\6» 1 be true for n ^ 1 so that the polar angle will stay 
well outside the range of logarithmic variation for high order modes. 

The distance at which a particluar higher order mode can be neglected is given by 



0a» a/\ v(n) \ — hjirn for n ^ 1 ((8) used for \v\ ). 



(12) 



Obviously the first order mode is the limiting case. The minimum distance is simply required 
to be large compared to the height of the ionosphere. (As in ordinary waveguides the modes 
below cutoff are negligible a few guide widths downstream.) A distance of 500 km satisfies this 
last requirement under either day or night conditions. Note that this result is independent of 
frequency (to the first order) over the range of 5 Hz to 1 kHz. 

An example of the relative magnitude of the zero order and the first order modes is shown in 
figure 1. The first two terms in the series for the normalized vertical field, E r l(r)Idsl2kha 2 ), are 
shown there. The zero order mode (n = 0) is shown in approximate form over the three ranges 
of distance that were outlined previously. The frequency of propagation is about 11 Hz, which 



Figure 1. Magnitude of the zero order and first order 
mode terms from the series for normalized vertical elec- 
tric field versus distance from the source, using approxi- 
mations from the text. 




FREQUENCY = 11 Hz 
IONOSPHERIC HEIGHT = 50 km 
SEE EQUATION (9) FOR DETAILS. 




L8= 11 
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corresponds to a minimum in the natural resonance structure. The height of the ionosphere is 
50 km as might be appropriate during the daytime. More will be said about the behavior of the 
zero order mode (with distance from the source) in later parts of this paper. But note that the 
first order mode has no important relative magnitude until the distance from the source is below 
200 km. The selection of 500 km as the minimum distance for predominance of the zero order 
mode is very conservative. 

3.2. Bounding Distances Between the Three Approximations for the Zero Order Mode 

Figure 1 demonstrates the point that for low frequencies there is a range of source distances, 
from 500 km out to the place where | v\ 6= 1, over which the normal asymptotic form for P v {— cos 0) 
does not apply. (The approximations of sec. 2 have been used in plotting the three ranges of dis- 
tance shown for the zero order mode.) The equivalent distance from the source that lies between 
these two ranges of approximation, is equal to 



0a = al\ HO) | = al(- (l/2)+V(A;aSo) 2 + (l/4)) 



(13) 



where use has been made of (5). Figure 2 is a plot of this boundary distance from the source as 
a function of propagation frequency. If the distance is measured instead from the antipode of 
the source then the same curve is a second boundary between a second pair of distance ranges 
based on | v \ (n — 0) = 1 (see fig. 1). Note that at higher frequencies this boundary is closely given 
by one-sixth wavelength. 

It should be emphasized that the necessity for three different types of approximation to 
P v {— cos 6) really only arises at the lower end of the ELF band in the normal propagation analysis. 
Figure 2 shows that the range of 5 Hz to 100 Hz is the frequency band affected. When the fre- 
quency exceeds 100 Hz the boundary distance is less than 500 km. The boundary distance is 
exceeded before the 500 km point is reached at which the zero order mode predominates. Thus, 
it can be ignored in the normal propagation analysis that assumes solely zero order mode. Strictly 
speaking the special approximation for the area around the antipode of the source still applies. 
However, one might argue that the area involved is small enough for frequencies above 100 Hz 
that it could be ignored for all practical purposes. 

Near the low end of the ELF band at 5 Hz the curve of figure 2 shows another phenomenon. 
The boundary distance between approximation regions has reached one-quarter of the distance 
around the earth. This boundary applies both around the source and its antipode. This means 
that the middle range of distances has disappeared, and the normal asymptotic expansion for 
Pv(— cos 6) does not apply at any distance. 

3.3. Higher Order Modes Within 500 km of the Antipode of the Source 
One point of rigor remains in the treatment presented here. The antipode of the source is 
a natural focus at which the various mode amplitudes can be expected to increase. The zero 



\ONE QUARTER CIRCUMFERENCE, 6 = 90 



MINIMUM DISTANCE AT WHICH THE 
ZERO ORDER MODE PREDOMINATES 



FREQUENCY IN Hz 



Figure 2. Distance from the source that separates two 
ranges of approximation for the zero order mode. 

(Same curve is the boundary between another pair of ranges where distance 
is from the antipode of the source.) 
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order mode should be checked to see that it predominates here also. Equation (10) shows the 
magnitude of the higher order modes would grow without bound except that that result does not 
apply within 500 km of the antipode. (An argument similar to that of sec. 3.1, but using 
| v(n) | (tt — 6) » 1 instead of \v{n)\d» 1, shows that (10) does not apply within 500 km of the 
antipode as well as 500 km of the source.) The correct approximation for high order mode terms 
within 500 km of the antipode is 

| v(y + l)P v (— cos 0)/sin V7T | — 2(iranlh) 2 exp (— (7ran/h)7r) 

for n ^ 1 

if | v | (77 — 6) « 1 (see sec. 2 for approximations). (14) 

The 6 dependence has disappeared as the mode amplitude approaches a constant. The net 
conclusion is that even around the antipode the higher order mode terms (n ^ 1) are negligible 
compared to the zero order mode. The exponential factor of (14) insures this fact just as in (10). 
The large differential attenuation between the higher order and zero order modes is not overcome 
by the focusing experienced at the source antipode. 



4. ELF Noise Spectrum Extensions 

The integral for ELF noise power spectrum from a vertical stub antenna, due to lightning dis- 
charge, has already been derived [Galejs, 1961a]. The section evaluates this integral expression 
when the sources are near the antenna or its antipode and the low-frequency end of the ELF band 
(5 to 100 Hz) is being considered. For the readers convenience several parameters are redefined 
in this section. A similar discussion, appropriate to the vertical electric field intensity has been 
given by Wait [1960b]. 

4.1. Integral Expression 

The starting point is the integral, just referred to, for received noise power spectrum of a short 
vertical antenna. 

rr \ ^- , 27n^+l)o-(o-+l) 

G(ico) = g(ia)) ,,,2 0,0. — t— — (15) 



1 



16k 2 a 2 h 2 sin per sin <tit 

P„(-cos 0)Pa{- cos 6) sin Odd 

0i 



where g(io)) = quantity proportional to the squared dipole moment of the lightning discharges per 
unit area (assumed independent of position), and 

co = radian frequency, 

k = free space, plane wave propagation factor, 
a = radius of the earth, 

h = height of the lower edge of the ionosphere, 

8 = polar angle at the center of the earth measured from noise source to receiving antenna, 
So = zero order root of the modal equation, a complex number that is largely real and close to 
unity. 

*^+l) = (A;aS ( )) 2 , 

<r = complex conjugate of v, 
P v {z) = Legendre function of z with degree v. 



767-936 0-65— 6 
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This expression assumes a constant surface density excitation of the earth-ionosphere cavity 
over a belt bounded by the two polar angles, 6\ and #2. Variations of surface excitation are ac- 
commodated by a piecemeal application of (15). Keep in mind that the radio-propagation depends 
only on polar angle and not azimuth angle. 

In this paper approximations will be made after integration so that all three possible ranges 
of approximation can be handled with one integration. According to Bateman Manuscript Proj- 
ect [1953] the integral equals 



/. 



02 

P v {~ cos d)P*{- cos 6) sin Odd 



= [" sin fl^-cos 0)Pl(-cos6)-PA-cos 6)P 1 fr (- cos 0)) J 2 
L {v-(j)(y + cr+\) J*i 



There does not seem to be any simple reduction of this form. So the next step is to approximate 
the various Legendre functions to obtain a result in convenient form for numerical computation. 



4.2. Approximate Results 
Consider first the case when the polar angle 6 lies in the middle range. Using the asymptotic 
forms of (2), combined with (15) and (16), the noise power spectrum is approximated by 

re \ «r ^ 1*^+1) I 2 

G(i(x)) = guar 



i(kah) 2 1 v + (1/2) I (cos (v - o)ir - cos (v + &)tt) 

cos {v + a + 1) (77 — 9) sin (v — cr)(7r— 0) 



02 



(^ + cr+l) (v-&) 

for \v\6» 1 and \v\{ir-0)» 1. (17) 

This is equivalent to the answer obtained previously by Galejs [1961a]. It has already been 
noted that the asymptotic forms used here for P v and P\ break down when the source is close to 
the antenna or its antipode. 

When the noise sources are near the antipode of the antenna then the approximations of (3) 
are used with (15) and (16) to get 

G(i<o)~g(i(o) An nT7 , *"t * / , x x -[-sinflcos 0/2]!! 2 for \v\(it-6)«1. (18) 

& ^{kahf (cos (^ — cr)7r — cos (v + a)rr) L J01 ' ,v { } 

The noise power observed drops to zero as the source approaches the antipode due in part to the 
vanishing circumference of the belt being integrated. (This is contrary to the middle range ex- 
pression that is finite at the antipode.) 

The remaining condition is when the sources are near the antenna, itself. (But still farther 
than 500 km away so that only the zero order mode is present.) Then the approximations of (4) 
are used with (15) and (16) to obtain 

G(ico) =« g(i(o) - 



2{kah)\v- &)(y + cr+1) 

. ( sin(^-o-)7T (v-or) ^ 1 1 r » m$2 

1 (cos (^-o-)7r-cos {V+&M it 2, | v + 1 + n\A ' LC ° S ' ^ 

(19) 

for \v\e<\. 
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FIGURE 3. Magnitude of ratio error squared from using 
midrange approximation, near the antenna or its anti- 
pode, for the electric field intensity. 




2 .4 6 14 16 18 20 

DISTANCE IN THOUSANDS OF KILOMETERS FROM THE SOURCE 



Thus, when the sources are near the antenna itself the contribution to the observed noise power 
spectrum still remains finite. This is again because the vanishing circumference of the belt being 
integrated tends to dominate the process. 

These three results for noise power spectrum [(17), (18), and (19)] cover all the possibilities. 
It may be that the limits of integration are in different regions of approximation. Since the exact 
integration is being approximated such a mixture of limits is easy to handle. Each limit is used in 
the appropriate approximation and the two subtracted to get the total answer. 

In figure 3 plots are presented to show the large error near the antenna or its antipode in cal- 
culating noise spectrum contributions using only the midrange approximation. The same param- 
eters are used in figure 1; frequency is about 11 Hz and ionospheric height equals 50 km. What is 
plotted in figure 3 is essentially the ratio of noise power contributions from an elementary sized 
area as approximated these two different ways. 



5. Discussion and Conclusions 

It has been shown that for terrestrial radio wave propagation in the middle ELF band (fre- 
quencies of 5 Hz to 1 kHz) the zero order waveguide mode is predominant when the distance from 
the source exceeds roughly 500 km. The theory of noise power spectrum in the ELF band that 
is extended in this paper assumes a "locally quiet" receiving site. In the sense employed here 
the phrase "locally quiet" means there are no thunderstorms within 500 km. Only the zero order 
modes propagating from lightning strokes need be accounted for. 

It has also been shown how the theoretical terrestrial ELF noise spectrum has different ap- 
proximate forms depending on whether the source is close to the antenna, or its antipode, or 
neither. This last modification of theory is chiefly important for frequencies below 100 Hz. The 
critical distances from the source or its antipode are roughly given by one-sixth wavelength. 

It is interesting to look at plots by Wait [1960b] of log (\E z \vp) versus (p) for frequencies down 
to 50 Hz. (See fig. 1 of that reference; E z is the vertical electric field and p is the distance from the 
noise source or lightning stroke.) At 50 Hz the curve is linear clear down to a distance from the 
source of less than 400 km. Presumably this is due to the predominance of the zero order mode at 
larger distances. This plot at 50 Hz shows that Wait's restriction of one-sixth wavelength for the 
limiting distance for linearity is more stringent than necessary. (At 50 Hz, a distance of one-sixth 
wavelength equals 1000 km.) 
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The zero electric field at 6= 112.5° in figure 1 is a result predicted by the same reasoning that 
accounts for the ELF resonance and antiresonance frequencies. As Wait [1964] has observed 
the resonance frequencies are closely given by 

(ft | So | ) (27ra) - 7T ~ 27TTI (20) 

where k = free space, plane wave propagation factor 
So = zero order root of the modal equation 
a = radius of the earth 
n = l, 2, 3. . . . 

The explanation is that resonance occurs when the magnitude of phase retardation for one circuit 
around the world is equal to an integer number of complete cycles. The subtraction of ir accounts 
for separate tt/2 phase advances that occur at the caustic points of each source and its antipode. 
Similarly the antiresonance points occur when 

(k\So\)(27ra)-7T^27Tn-7T (21) 

so that a complete circuit tends toward phase cancellation. 

Following the same train of thought there should be minimum electric field intensity at an 
angle from the source 6 where 

[(ft | So | )((2tt - 6)a) - (3tt/4)] - [(kS ) (0a) - (tt/4)] - 7r(2n + 1). (22) 

The difference in phase retardations equals a number of complete cycles plus one-half. In the 
example plotted in figure 1 the frequency of 11 Hz is an antiresonant point where ft|S |a — 2. 
Substituting this in (22) and solving gives 9 — 112.5° as the only angle in the middle range where 
phase cancellation tends to occur. The graph verifies this conclusion. 

The extensions of the formulas for terrestrial noise power spectrum given in this paper will 
have application in a future paper concerned with accounting for temporal and geographical 
distribution of lightning strokes. 
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